AD-A078  655 
UNCLASSIFIED 

To?8665  I 


INDIANA  UNI V  AT  BL00MIN6T0N  DEPT  OF  CHEMISTRY  R/6  7/4 

AN  ACCURATE  MODEL  FOR  SAMPLE  DROPLET  ACCELERATION  IN  ANALYTICAL— EtC (U) 
DEC  79  R  E  RUSSO  »  G  M  HIEFTJE  N00014-76-C-0838 

TR-24  NL 


END 

DATE 

FILMED 

I  -80 


DOC 


Security  classification  or  this  p*cf  '»/, 


RKAU  INSTRUCTIONS 
BEFOK'i:  COMPLETING  FORM 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  number 


ENT'S  CATALOG  NUMUEN 


2  GOVT  ACCESS! 


TWENTY -FOUR 


An  Accurate  Model  for  Sample  Droplet  Accelerati 
in  Analytical  Flames  ,  »  - - — - - 


Interim 


».  AUTHORfAj 


•  .  CONTRACT  OR  GRANT  NUMBER/*,) 

N-14-76-C-0838  / 


(0*  Richard  E./ Russo 


PERFORMING  organization  name  ANO 
Department  of  Chemistry  / 
Indiana  University 
Bloomington,  IN  4740S 


WfcGRAM  ELEMENT.  PROJECT 
REA  A  WORK  UNIT  NUMBERS 


CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Office  of  Naval  Research 
Washington,  D.C. 


MONITORING  AGENCY  NAME  ft  ADDRESS///  d///«j 


IS.  SECURITY  CLASS,  /o/  //■/•  r«port) 


UNCLASSIFIED 


1S«,  DECL  ASSIFICATION/OOWNGRAOING 

schedule 


i*.  distribution  state“nt  to/  »m,  Report) 


Approved  for  public  release;  distribution  unlimited 


IMENT  (ofphe  Jhat&et  entered  In  Block  20,  II  different  from  Report) 


1§.  SUPPLEMENTARY  notes 


Prepared  for  publication  in  ANALYTICAL  CHEMISTRY 


KEY  WORDS  (Continue  on  reverse  tide  If  necessary  and  Identity  by  block  number) 


24.  ABSTRACT  (Continue  on  reverse  sld e  It  neceseery  end  Identify  by  block  number) 


The  experimentally  observed  velocity  of  an  aerosol  droplet  in  an  aii -acetylene 
flame  is  compared  to  the  velocity  calculated  through  a  new  mathematical  model. 
The  experimental  droplet  velocity  is  measured  by  means  of  stroboscopic  p  y  >g- 
raphy  of  a  stream  of  individual  droplets  sent  into  a  flame  at  a  known  frequency 
and  traveling  very  reproducible  trajectories.  Because  the  temporal  spacing 
between  these  droplets  is  known,  a  measure  of  the  instantaneous  droplet  i  ca¬ 
tion  enables  their  velocity  to  be  calculated.  The  droplet  acceleration  ,  >de 1 
was,  modified,  bv  incomorat inc  into  it  a  snatially.  changing,  flame  ris" - - - - : 


form 
1  JAN  73 


EDITION  OP  «  NOV  65  |S  OBSOLETE 
S/N  0  102-014-  6601  y 


UNCLASSIFIED _ 

TV  CLASSIFICATION  q  r  TMiS  PACE  rwi.rn  Ds'e  Enter**) 


UNCLASSIFIED 


tUlHIty  CUHiriCUlON  OF  THIS  P»c,ff*li«i  Dmim  Kn"t,d) 


20.  (Abstract  -  continued) 


"’/Velocity,  which  was  assumed  constant  in  earlier  treatments.  Excellent 
agreement  has  been  found  between  calculated  and  measured  droplet  ac¬ 
celeration  in  the  flame.  k 


StCUWTY  CLASSIFICATION  OF  THIS  P  AGI.fWTiOf  C«r»  I  > 


OFFICE.  OF  NAVAL  RF.SEARCU 
Contract  N14-76-C-0838 
Task  No.  NR  051-622 
TECHNICAL  REPORT  NO.  24 

AN  ACCURATE  MODEL  FOR  SAMPLE  DROPLET  ACCELERATION 
IN  ANALYTICAL  FLAMES 


by 

Richard  L.  Russo  and  Gary  M.  Flieftj 


Prepared  for  Publication 
in 

ANALYTICAL  CHEMISTRY 

Indiana  University 
Department  of  Chemistry 
Bloomington,  Indiana  47405 

December  5,  1979 


\  Accession  F°r  a 

\  NTIS  fiRAfcl 

; 

1  Un  a-  •  a  " 

l  - 

Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government 


Approved  for  Public  Release;  Distribution  Unlimited 


BRJfcF 

wow 

A  revised  model  for  droplet  acceleration  in  an  analytical  flame  or 
plasma  accurately  predicts  experimentally  observed  behavior  in  an  air- 
acetylene  flame.  The  principal  modification  involves  a  spatially  changing 
flame  rise  velocity. 


ABSTRACT 

vwwwv. 

The  experimentally  observed  velocity  of  an  aerosol  droplet  in  an  air- 
acetylene  flame  is  compared  to  the  velocity  calculated  through  a  new  mathe¬ 
matical  model.  The  experimental  droplet  velocity  is  measured  by  means  of 
stroboscopic  photography  of  a  stream  of  individual  droplets  sent  into  a 
flame  at  a  known  frequency  and  traveling  very  reproducible  trajectories. 
Because  the  temporal  spacing  between  these  droplets  is  known,  a  measure 
of  the  instantaneous  droplet  location  enables  their  velocity  to  be  cal¬ 
culated.  The  droplet  acceleration  model  was  modified  by  incorporating 
into  it  a  spatially  changing  flame  rise  velocity,  which  was  assumed  con¬ 
stant  in  earlier  treatments.  Excellent  agreement  has  been  found  between 
calculated  and  measured  droplet  acceleration  in  the  flame. 
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INTRODUCTION 

WjVWWWI, 

An  aerosol  droplet  of  sample  solution,  directed  into  an  analytical 
flame  or  plasma,  will  rapidly  accelerate  toward  the  moving  gas  velocity. 

The  rate  and  degree  of  approach  of  the  droplet  to  the  flame  gas  velocity 
is  important  both  fundamentally  and  practically  in  analytical  spectrometry. 
For  example,  the  differential  velocity  between  the  droplet  and  the  flame 
(or  resulting  particle  and  the  flame)  governs  whether  or  not  convective 
mass  transport  is  significant  in  desolvation  and  vaporization.  Also,  drop¬ 
let  acceleration  affects  the  eventual  distribution  of  analyte  atoms  through¬ 
out  the  flame,  which  in  turn  influences  analytical  sensitivity,  precision 
and  working  curve  linearity. 

In  particular,  the  horizontal  acceleration  of  droplets  as  they  enter 
a  flame  has  been  shown  to  be  responsible  for  some  of  the  lateral  spread 
of  atom  populations  (l).  More  importantly,  the  vertical  acceleration  of 
aerosol  droplets  determines  the  distribution  of  analyte  atoms  above  the 
burner,  and  therefore  contributes  to  the  optimal  viewing  height  in  the 
flame.  It  is  also  this  vertical  acceleration  which  determines  the  amount 
of  time  before  a  droplet  (or  particle)  reaches  the  flame  gas  velocity, 
at  which  time  convective  mass  transport  becomes  negligible  in  its  influence 
on  the  atomization  processes. 

A  mathematical  expression  has  been  formulated  to  model  the  vertical 
acceleration  of  aerosol  droplets  in  a  flame  or  plasma  (2).  However,  the 
time-dependent  droplet  velocity  calculated  through  this  expression  docs 
not  accurately  predict  the  velocity  actually  observed.  The 


apparent  oversight  in  this  previous  work  (2,3)  is  that  fl ante  rise  velocity 


was  assumed  to  be  constant  across  the  burner.  Unfortunately,  this  assumption 


leads  to  greater  theoretical  than  observed  droplet  acceleration;  in  the  ex¬ 
perimental  test  ot'  the  mode  1,  droplets  were  injected  into  the  flame  in  a 
direction  perpendicular  to  its  flowing  axis,  so  they  experienced  a  spatially 


changing  flame  velocity.  Moreover,  in  the  present  study  it  was  found  that 


the  measured  droplet  velocity  in  the  flame  is  affected  significantly  by 


the  droplet  introduction  rate.  At  droplet  introduction  rates  greater  than 


2KHz,  the  droplets  bunch  together,  forming  a  quasi -cont inuous  stream  of 


solution.  The  flame  then  must  accelerate  the  entire  stream  instead  of 


individual  entities,  therein  lending  to  different  behavior  than  expected 


for  acceleration  of  s 


In  the  present  paper,  the  agreement 


between  theoretical  and  experimental  results  is  improved  by  incorporating 


a  measured,  spatially  changing  flame  rise  velocity  into  the  acceleration 


expression,  and  by  controlling  carefully  the  droplet  introduction  rate 


Significantly,  theoretical  results  are  shown  to  compare  quite  favorably 


with  the  measured  droplet  velocity  over  the  entire  droplet  lifetime  in  the 


TUTOR Y 
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The  droplet  acceleration  expression  is  derived  from  Newton’s  Second 


Law,  a=F/m,  where  a  is  the  acceleration  of  the  droplets,  m  ; s  their  mass,  and 


F  is  the  algebraic  sum  of  forces  acting  on  the  droplets  in  the  flame.  There 


& 


to  accelerate  the  droplets,  and  gravitational  attraction,  which  tends  to 
retard  them.  For  a  laminar  flame  and  spherical  droplets,  viscous  drag 
can  be  expressed  by  Stoke's  Law  as:  3*nd(V  V,),  where  n  is  the  average 
viscosity  of  the  flame  gases,  d  is  the  diameter  of  the  droplets,  and  V 
and  Vj  are  the  velocities  of  the  flame  and  droplet,  respectively.  Grav¬ 
itational  attraction  is  mg.  where  g  is  the  gravitational  constant. 

Unfortunately,  droplet  diameter  and  mass  are  not  constant  in  the  flame 
because  of  ongoing  desoivation.  Therefore,  time-dependent  relations  for 
droplet  mass  and  diameter  must  be  incorporated  into  the  foregoing  treat¬ 
ment.  The  time  dependence  of  droplet  sire  during  desolvation  has  been 
determined  (4)  to  be,  d1  -  d^  -  kt .  where  d  is  the  droplet  diameter  at  time 
t  after  its  entry  into  the  flame,  d^  is  the  initial  droplet  diameter,  and 
k  is  the  desolvation  rate  constant.  Upon  substitution  of  the  two  time- 
dependent  forces  into  Newton's  Second  Law,  the  droplet  acceleration  ex¬ 
pression  becomes: 


_  isn(v  -  va) 

'p(dz  -  kt 


g 


where  p  is  the  density  of  the  aerosol  droplet. 

In  an  earlier  treatment  (2),  this  expression  was  integrated  to  calculate 
the  droplet  velocity  in  the  flame.  However,  direct  integration  is  valid  only 
if  one  assumes  the  flame  velocity  to  be  constant  over  the  entire  lifetime  c^" 
the  droplet.  We  have  found  that  this  assumption  is  not  valid,  especially 
when  a  droplet  generator  (41  is  used  to  iniect  the  aerosol  droplets  into 
the  flame.  In  this  situation,  used  to  tost  the  proposed  model,  droplets 
are  introduced  into  the  flame  in  a  direction  perpendicular  to  the  flame 
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flow  and  therefore  experience  a  parobolic  flame  rise  velocity  gradient  as 
they  penetrate  into  the  flame.  Consequently,  to  obtain  agreement  between 
the  proposed  model  and  results  obtained  using  the  droplet  generator,  one 
should  determine  experimentally  the  flame  velocity  profile,  relate  this 
profile  to  the  time-dependent  horizontal  droplet  position,  and  then  inte¬ 
grate  the  acceleration  expression.  Unfortunately,  such  a  procedure  leads 
to  a  fourth-order  dependence  of  acceleration  on  flame  rise  velocity,  making 
integration  difficult  and  laborious.  Instead,  a  numerical  approach  was 
chosen  in  the  present  study.  In  this  approach,  the  acceleration  expression 
is  directly  evaluated  at  very  small  time  intervals  using  an  experimentally 
mapped  flame  rise  velocity. 


EXPERIMENTAL 

VWWWAWW, 

In  the  measurement  of  droplet  velocities,  a  droplet  generator  is  used 
to  inject  into  the  flame  uniform-sized  aerosol  droplets  with  reproducible, 
isochronal  spacing  (4,5),  Because  of  this  reproducibility,  a  single  high- 
intensity  flash  photograph  reveals  the  trajectory  and  spatial  relationsh  p 
of  droplets  in  the  flame,  from  which  their  instantaneous  and  spatially 
dependent  velocities  can  be  deduced.  Figure  1  of  reference  4  shows  such 
a  single-flash  photograph,  with  a  droplet  introduction  rate  of  1000  drop¬ 
lets  per  second. 

In  the  present  study,  droplets  were  introduced  into  the  flame  in  a 
slightly  downward  direction  in  order  to  establish  a  zero-vcrtical-velocit 
droplet  location.  Experimentally,  vertical  droplet  positions  were  o''tr  .nod 


i  .■  ■ 
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by  direct  measurements  above  this  reference  droplet.  Of  course,  small 
variations  in  individual  droplet  positions  exist  because  of  air  currents 
about  the  apparatus.  Therefore,  droplet  positions  from  several  photographs 
were  averaged  to  obtain  an  experimental  curve  of  droplet  position  versus 
droplet  residence  time  in  the  flame  (Figure  1).  Each  droplet  position  in 
Figure  1  represents  the  average  of  five  different  photographs  in  which  all 
droplet  generator  and  flame  conditions  were  the  same.  Error  bars  in  Figure 
1  represent  95%  confidence  limits,  indicating  that  less  than  6%  error  exists 
in  assigning  droplet  positions.  If  necessary,  the  magnitude  of  this  error 
could  be  decreased  by  averaging  additional  photographs  or  by  shielding  the 
experimental  apparatus  from  ambient  air  currents. 

The  time-dependent  horizontal  position  of  th<  droplets  in  the  flame 
was  also  measured  from  the  single- flash  photographs.  Coupled  with  knowledge 
of  the  rise  velocity  of  the  flame  as  a  function  of  position  across  the 
burner  (horizontal  flame  rise  velocity  profile),  this  measurement  enables 
one  to  calculate  the  flame  velocity  that  each  droplet  experiences  as  it 
traverses  the  flame.  The  method  of  obtaining  flame  rise-velocity  profiles 
is  described  elsewhere  (6).  The  rise-velocity  profile  for  the  air-acetylene 
flame  used  in  this  study  is  shown  in  Figure  2.  Asterisks  in  Figure  2 
represent  the  flame  rise  velocity  values  measured  at  a  particular  burner 
location.  Again,  error  bars  show  the  95%  confidence  limits.  Not  suiprisingly , 
the  relationship  between  rise  velocity  and  position  across  the  flame  was 
found  to  fit  a  second-order  polynomial  expression  with  a  correlation  co¬ 
efficient  of  0.98.  The  smooth  parabolic  curve  (dotted),  obtained  from 
this  polynomial  expression,  is  what  one  would  expect  from  Poiseuille's 
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equation  for  a  laminar  flow  system.  The  velocity  of  the  flame  that  each 
droplet  experiences  with  respect  to  spatial  position  can  he  easily  found 
from  this  polynomial  (or  from  Fig.  2),  and  direct  evaluation  of  the  accel¬ 
eration  expression  at  short  time  intervals  then  becomes  possible. 

A  lOOppm  barium  chloride  solution  (density  =  1.0  g/inL)  was  used  in 
the  droplet  generator  reservoir.  The  initial  droplet  diameter  was  measured 
to  be  68  t  2  pm,  using  the  MgO  technique  (7).  Droplet  introduction  rates 
were  kept  below  2KHz  in  order  to  minimize  the  slip-stream  effect  (8).  Flame 
gas  viscosity  was  calculated  to  be  660  micropoise,  using  statistical  mech¬ 
anics,  for  a  flame  temperature  of  2450K  (9).  This  average  flame  tempera¬ 
ture  was  measured  over  the  flame  region  of  study  using  the  sodium  line- 
reversal  technique  (10). 

RESULTS  AND  DISCUSSION 

WVVVVWWVlA/WWWWV 

The  random  spatial  distribution  of  droplets  produced  by  a  conventional 
nebulizer  makes  it  unlikely  that  one  droplet  will  strongly  affect  the  ac¬ 
celeration  or  velocity  of  another.  However,  with  the  droplet  generator  used 
in  the  present  study,  droplets  follow  identical  trajectories  through  the 
flame  and  can  affect  each  other  by  developing  a  slip-stream  (8).  Under¬ 
standably,  this  slip-stream  effect  is  most  marked  when  droplets  are  closely 
spaced;  visual  changes  in  droplet  velocity  (spacing)  are  evident  at  dropl  t 
introduction  frequencies  greater  than  2Kl!z.  Because  this  slip-stream  ef¬ 
fect,  if  pronounced,  tends  to  shield  each  droplet  somewhat  from  the  viscous 
forces  of  the  flowing  flame  gases,  upward  droplet  acceleration  is  somewha* 
lowered  at  high  droplet  introduction  rates. 


This  reduction  in  acceleration  is  evident  by  examination  of  Figure  3, 
in  which  all  experimental  conditions  were  the  same  except  for  the  droplet 
introduction  rate.  At  an  introduction  rate  of  2.8KIIZ,  droplets  approach 
the  flame  velocit>  more  slowly  than  at  l.GKHz.  To  obtain  a  measure  of 
droplet  velocity  consistent  with  the  developed  theoretical  model,  a  low 
droplet  introduction  rate  (l.GKHz)  was  employed  in  the  present  study.  At 
this  rate,  the  slip-stream  effect  should  be  negligible  (8). 

Figure  4  compares  droplet  velocity  curves  calculated  through  the  ac¬ 
celeration  expression  with  one  obtained  experimentally  at  a  droplet  intro¬ 
duction  rate  of  l.GKHz.  The  upper  theoretical  curve  (dot-dash)  represents 
the  droplet  velocity  calculated  assuming  a  constant  flame  rise  velocity  of 
10.0  m/s.  The  lower  curve  (dotted")  was  calculated  by  incorporating  the 
spatially  changing  flame  rise  velocity  (ca.  Fig.  2)  into  the  acceleration 
expression  and  evaluating  it  as  small,  discrete  time  intervals.  As  expected, 
the  predicted  acceleration  is  too  great  when  a  constant  average  flame  rise 
velocity  is  assumed.  In  contrast,  the  curve  incorporating  a  spatially 
changing  flame  rise  velocity  (dotted)  agrees  within  experimental  error 
with  observed  velocities  over  the  entire  lifetime  of  a  desolvating  drop¬ 
let.  Clearly,  a  spatially  changing  flame  rise  velocity  must  be  considered 
in  order  to  accurately  model  droplet  acceleration  when  a  droplet  generator 
is  used  for  nebulization.  In  conventional  nebulizer  systems,  however, 
droplets  usually  enter  the  flame  vertically  and  a  constant  flame  rise 
velocity  in  the  acceleration  expression  should  provide  close  agreement 
with  observed  behavior. 

Additional  but  small  imporvement  in  the  droplet  acceleration  model  is 
possible.  For  example,  the  existing  expression  does  not  take  into  accou-'t 
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the  change  in  density  which  occurs  as  droplets  desolvate  and  become  more 
concentrated.  Also,  a  changing  flame  viscosity  could  be  incorporated 
into  the  model.  Such  a  change  probably  occurs  as  a  result  of  changing 
flame  temperature  and  water  vapor  concentration  in  the  vicinity  of  the 
drop  ts.  Inspection  of  the  acceleration  expression  reveals  that  both 
of  these  modifications  would  improve  agreement  with  experiment.  Finally, 
o.t  is  not  yet  certain  that  the  effect  of  one  droplet  on  another  is  completely 
removed  at  the  droplet  introduction  rate  employed  in  this  work.  In  order 
to  measure  the  true  droplet  trajectory  (velocity)  in  the  flame,  one  would 
have  to  inject  a  single  droplet  and  monitor  its  path  with  a  high-speed 
motion  camera. 

In  our  view,  the  improved  agreement  between  experiment  and  theory 
which  the  above  three  modifications  would  yield  does  not  justify  the  sub¬ 
stantial  complications  which  they  would  introduce  into  both  the  existing 
model  and  its  verification.  In  most  cases,  the  agreement  obtained  in  this 
study  would  be  sufficient  to  predict  both  the  spatial  distribution  of  atoms 
in  an  analytical  flame  and  the  amount  of  time  required  for  introduced  drop¬ 
lets  (or  particles)  to  reach  the  flame  rise  velocity. 

CRliPri 
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Figure  1.  Droplet  trajectory  in  the  flame.  Initial  droplet  diameter 
*  68  um.  Introduction  rate  -  1.6MU.  Error  bars  indicate 
i  95%  confidence  limits. 

Figure  2.  Spatial  velocity  profile  in  stoichiometric  (fuel/oxidant  *  0.15) 
air/acetylene  flume.  Asterisks  represent  measured  velocities; 
dotted  line  is  curve  obtained  from  second-order  fit  (correlation 
eofficient  =  0.98). 

Figure  3.  Measured  droplet  velocity  in  the  air  acetylene  flame.  All 
droplet  generator  and  flame  conditions  are  the  same,  except 
for  the  droplet  introduction  rate. 

0  introduction  rate  =  2.8Kllz 

*  introduction  rate  »  l.hkllz 

Figure  4.  Comparison  of  measured  droplet  velocity  with  that  calculated 
through  the  acceleration  expression  (hq.  1). 

*  =  measured 

calculated  assuming  rise  velocity  is  constant  across 
flame. 

-  calculated  using  horizontally  changing  flame  rise 
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